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svstf.m- aivn m fthod for analyzing and imaging 
Three-dimensional Volume DataS^is 



Background of the Invention 



Field of the Invention 

The present invention relates generally to imaging of three-dimensional ("3D") 
volume data sets. More particulai'ly , the present invention relates to fast visualization and 
analysis of structures within 3D volume data sets. 

Related Art 

Many fields of endeavor require the analysis and imaging of three-dimensional 
("3D") volvime data sets. For example, in the medical field, a CAT (computerized axial 
tomography) scanner or a magnetic resonance imaging (MRI) device is used to produce 
a "picture" or diagnostic image of some part of a patient's body. The scanner or MRI 
device generates a 3D volume data set that needs to be imaged or displayed so that 
medical persormel can analyze the image and form a diagnosis. 

Thi'ee-dimensional volume data sets are also used in various fields of endeavor 
relating to the earth sciences. Seismic sounding is one method for exploring the 
subsmface geology of the earth. Animdergroimd explosion or earthqualce excites seismic 
waves, similar* to low frequency sound waves, that travel below the surface of eaith and 
are detected by seismographs. The seismogi-aphs record the time of arrival of the seismic 
waves, both direct and reflected waves. Knowing the time and place of the explosion or 
earthquake, the time of travel of tlie waves through the interior can be calculated and used 
to measure the velocity of the waves in the interior. A similar- technique can be used for 
offshore oil and gas exploration. In offshore exploration, a ship tows a sound soxirce and 
xmderwater hydrophones. Low frequency (e.g., 50 Hz) sound waves are generated by, 
for example, a pneumatic device that works like a balloon burst. The sotmds bounce off 
rock layers below the sea floor and are picked up by the hydrophones. In this manner, 
subsurface sedimentary structures that trap oil, such as faults, folds, and domes, are 
"mapped" by the reflected waves. The data is processed to produce 3D volume data sets 
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, f„, some measured or calculated property. e.g..seisnno amplitude 
.asanu^ene,™ - ^^^^^^^^ 

otthe volume atthatlocatton. O plurality ottwo- 

3D volume data set IS to cro.s-sert,onfl.e 3D ^^,^^,,,,3,Hstiiea 

• „»l ^"7n"^ cross-sections or slices. The image of tie 3D volume 
dmiensional ( 2D ) cross offlie 3D volume data set 

rrrce:^::rt~:aud.eucomposited.to 

'^--^'".^^*t:rrrl.e.eius«tlayer.y layer asthecomposite 

""'"^ ;,4,^,Lorga.u.ationorstruc^of*evo>umeas 

""rilslrcLdoompositeteclrni^ueistypicaUy 

r:::::r::...vo..ed.se..^^^^ — - 

teau^res i.ie«ut it, the 3D volume data set 

romouter software has been developed specifically for .magmg 
Computer sottw „„,,,of ^uoh conventional computer programs 

sets for the oil and gas industry. Examples of su^h 
^oludeVo.elOeo.avai,ahle.omParadismOeophystca^^— 

Bar^Cuhe. available .om Landmarlt Graphics ^^'^^^^^^r.^^^^ 

Tarriiciset. 

interpreting 3D seismic data operate o„ «'\*^ ^ „,ey or 

consequently, every thne a change is made, such as "-^^ssed and *e image 

.u. M\ volume of seismic data must be processeo, <uiu 
opacity settmgs, the Ml 3D volume 0 ^^^^^ 

0.02 frames per second, respectively). Durmgm 
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. on ^ of *res.. maHng it difficult to ccnpletely and 

the user loses focus on the teauii 

P"P='>^^'™*'rt;ol^cinte,«ta.onptogt^ 

visualize and in...p»t a piece offhe HI 3D vo ^i„age of fl.e selected 

elected uieoe via a tnenucoiranana an e 
coordinates of the selected pt ^ 

,.e is dra^n. The selected piece o. .» -^^^^ 

Ho«eveMoloo.a.adi«e«n.P-"J-^^^ „„3t .e 

«lo. a seoiosic feature ^^^^^ .^nt piece ts entered, and a new 

„ed. a new location or coordr^e ^ ^^^^^^ 

,„aae is drawn containing *e d.ffer». p.^^ ^ ^„ pieces of *e 

^es it difScnl. for tire user to vrsnaltze any _^ ^^^^ 

-SOvolnrneofseisrnicda...— ^^^^ 

interpret and identify the geologtcfeanare^tha „,^„f ^i^mic 

.isrracdat. Addi.iona,l.even.ho.ho.^ya^«^^-^^^^ 

....eingvisi^lydisplayeico.— -^.Jinr^ 
processing the full 3D votame of seisnnc data 

^'=^'-°^*°TCl^cin.erpretationprogra»spro^de*ecapa«lityto"anto 

'° r L r sa.s^ a voxel selection algotithnr. However. «tese 
pick" and identrft pomts fta. san ft- .^^^^^^ ^ 

programs t^ically iterate ^ough the Ml » v^ ^ 
p...^.satis.y.hevo.e,— .^^^^^^^ 

.peed graphics supercomputer. ^ — delete 6o. the coliecHon of picked 
p^o^ams do notprovide^e -P**-^^^ of picked voxels using 

voxels. Tl^eonly way to-elhninate" ^ ^lec.on 

^nventional 3D seisnric interpretation progran. '^ "J^^J^^^,^^,,^, 

...riaforthevoxelselectionalgorithnrunta^eP^--:^^^^^ 
..^electioncriteriaforthedisplayedp™..^-^^^^^^ 

^.htinre^eseiec^oncriterrats^^u-^^e^^^^^^^^^^^^^^^^ 
process is tinreconsnnnng, and — -^^^ ^^^^^^^ 
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WOOW^'W fwae 3D volume daasesm response to 

Ijin*e^»-«>*°"'*"^^* !r.e.^e^>ay.d image inaoon^nuo.^ 

.a, ..o« a use. .0 — * Sue. a ..e. 

l.e.>.ou...™P--^:l;.,.^ 

^ouldallowa.sertoxnorequicW> 
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,l3neandanopposln„rohef^P>-0^;-^^ ,^„»lpoln.^ 
leontrolpolntslntheprohefa^pla^-^ 

o^eorrno^lhtesondreprohet^ep^^- ^e,^^^^ 
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^^ensional sampling probe. to plurality of control pomts 
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„ more of pl^^l-tV 

locationofoneormoreofthepiur y 

^ , p^erred emboto-^then^b ^^^^^^^^^^^^^^^^^ 

.0 Ota plurality of side of ^^^^^ , „.eWng a particular 

I„3r,o*erembod»entof.hepre.e^^^ method ntay comprise *e steps of 

physical phenomena, sueh as a geo ogre . ^^^„,^i„„al volume 
p^itiomng^eprobef^P-— ^ 

^^..aadform-. fir^^s^of ^ 
15 pHysl-lP"— ^«=^"^"^*;^^^„t^,.„ipointstodefi«eafirs.spliue 
^, .elude interpolating 

or^e Other steps may meludemovmg .hep „^ ^re pmbe face 

.^e.ionalvolumedataset,formingas^nds.^^^^^^ 

pianefortrac^ngthephysicalphenomenaandm^^^o^ « ^^^^^^^ 

interpolating a three dmensional surface betwe 
,,necurve*hisrepresen««veof*ep«— 

-'-«^-*-^«-""r~md set of control points. Uisarr 
,..rsectsthefirsts.ofcon,ro,^^^^- 
25 advantageoffhepresentmventronttaltbeMs 

subs^uent spline curve, are curvilme.^ ^^^^ ^ 
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i spUne and the thkd spline emve tor turthei 
be«,een the first spline physical phenomena, 

defining the three dimensional surface ^ ^^^^ 

. des^e, the -X;::::^:eposifionsoftheprohe. 
.^presentWecontrolpo^tsont^e^oh^.^ ^^^^^ ^.^ ^_„, 

S Mor^ver. the me*od nr. ^^'X4„,,„^„,^.p„,.. Ano.ers.epma, 

^.weenr^pectivccontrolpomtsatrespe^ P 

..^edUplayinsthesplinecn^es-^"^ 
T^esplinec^esandthe.^^^- 
rhephysicalphenomen. T^^^dmc 

r:rrZ:rpl.e,.here.reshap..es.aceand.d.^^^^ 

.,e current spline enrve and the prior ^^^^„r.om.,^^^^^ 
Preferablythemethodmayalsomcludestepss 

,f hrtersecaons to *erehy reposi«on .he prohe face ^^^^^^ Lo. Points . 
, in^rsection The method also comprises sele^S JJ^*-- ^^^^ 
,^.y.posifion.eprohefacep^.^^^^^^^^^ 

Stated anoAer way, an embodiment or m ^.uin the three-dimensional 

posiUoning*eprohefacep.anea.ap.,^tt.ofposr«— 

— -^-^°rr:::nri:::: sphne c^r.. 

P'- " Xeandinterpolalinsh^^^ 

FEATUBESA^^D ADVANTAGES ,^.,,,,^,^00 section through a 3D sampling 

Itisafeatureofthepresent— ^ n^^^ 

p.hecanbec.eated.ed.awn,edited,^d — 
apluxalityoflinesthataxethenprojectedlhroughtheSDsa^plmgJ^ 

J drawn at angles offset ^o. tl.e coordinate system, such as an . y, z 
coordinate system, ofthe 3D sampUng probe. ,,^es in a 3D data volume 

T. • another feature of the present invention that structures m a iu 
30 It xs another feature .„i.eauickly mapped by selecting pomts 

set such as for instance geological structures, can be quickly mapp 
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• ti..^Dsamplmgprobe,whichpomtsmaythenbe 
ofinterestatapluralityoflocations«.theDsa.p^^^^^^ 

, 1 a volume data set, tne ^ 
.30s™plin.P»l==*^ousha3D^^^ ^^^^^^^^^^ 

^.toe «ft movement of *e 3D p„be is 

— ^ 3D sampUng P'"^/^™^; ^ , ^ fly so « the user 

^0 volume undated »ith va„tng deg^s^ P ^^^^ ^^^^ 

Uisatetherfea^eofthepresentu ^^^^ ^^^^ 

.al-tinte Wiethe change .n shape ot of * 

of the 3D sampltog prohe ^ ^ interaetively 

"''^^"*"\'^"^elthe^aoesofthe3Dsa»p>h.gpt^^ 
„tatea3Dsamplingprob=sothattheB^ „ein«« changing in leal-ttae with 
,.-^™..on,hefiy.soth..ens.^^ - J^^^^^^^ 

station of tte 3D samplmg ptobc. "j ,^^ ^4 varying degrees 

25 of trmspaienoy ontneuj 

,Ath«heiotationofthe3DsampHngprobe. JD sampling ptobe 

,Usyeta*^erfeah^of*epresent— ^^^^ 

^ ^ created arrd manipniated by the nser to dnecUy delete 
points that fall vrithin a certain datavatae range. 
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+-^r. ihi^ii a user can manipulate a 3D 

image a feature. interactively 
It is a further advantage of the present invention that a user 

, „f 3D seismic data. The present invention can belied to 
interpret iaxgevolmnes of 3D setsnuc ^^^^,^^,^«g.„usly 
,^eay3ndaeonra.e..iden.ifyd..tass.tes. ^^~v^^^^^^ ^^^^^^^ 

be used to sharply reduce 3D seisnuc project cycle Umes, 
,5 existirrgv-ells, and to locate additional reserves. 

Thepresentinventtonrsdescnbed^tth^ 

. • c MVe reference numbers mdicate identical or 
In the dra^vmgs, like referenc .^^^^^^^ ^^^^ 

20 elements. Additionally, the left-most digxt(s) of a 

..awinBinwhich^ereferencenumberfn.app^s^^^^^^^ ^^^^^^ 
FIG.l shows one embodiment of a software or program 

implementing the present invention; 

probe program ofthe present invention; 

25 . a function of datavalue; 

FIG. 3 shows a curve illustrating opacity as a functio 

;n„«tratine one embodunent tor 
FIG. 4 shows a flow diagram illustratmg on 

implementing the present invention; 
FIG.5 showsaflowdiagramofoneembodimentforchangmgadefault 
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J- ^ nf one embodiment for creating 
FIG. 6 shows a flow diagram of one em 

additional probes; „,„rnbe- 

in 3D space; . 
no. n s.o^ a flow dia^ o( one ' fo. ™.a« a probe 

While fixed in space; 

picking or seed selection; 
,XG. n lows one e».oto». of a c„n.pu.e. sys.e« sunable fo, «3e 

with the present invention; 
,.,3 :lan-n..e— ofaco„pu.e,svs.e.s..Me 

for use with the present invention; 

for use With the present invention; 

two of the probes inteiseoting each other; 
^O.U mustrates^eeprobesof^epresent — --probe. 

ttansparentontprobcandavolunae-renderedprobe. 
,,0.n uJ-aribbonsectiona^rdht^o^epresent — . 

.,0 « — s'^^el din^enstona. ^d and «.ee dintenstona, 
leerepresen.3«veofaph..a.pheno.enades»^^Va 

3.D volunte data set according to the present — , a.d 
^a. 20 shows a bloc, diagram of a systetn for producing *e grrd and 
surface of F1G.19. 
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I,..An.EI,D.SCKlPTIONO.THEPKEFEKKEI>EMBODIMENTS 

. mn is directed to a system and method for analyzing and 
The present mvention is directed 

• ai cnD"^ volume data sets using a sampimg p 
imaging three-dimensional ( 3D ) volum ^.^^oxel is a sample or point 

• 'S/nvpl <;" or volume elements. Eacn voxei is a ^ y 
volmnedatasetscompnse voxels orvoi datavalue) where "x, 

.ithinavolume. Eachvoxelcanheexpressedm^e orm^ M^^^^^^^ 

,.i.enti.esthe3mooationonhepoint^^-^^^ 
..ueofsomemeasuredorcalculate^a..^^^^^^ 

^^^^^^ 

waves at the specified (x, y, z) location. 

Any form of infomiation that can be expressed m the 
Any luiiii . x- „ Trt addition to seismic data, 

examples from the oil ana g .^^^^ 
^^^rip\t\t remote sensing data, ana siucav. 

— — -t"— ..example... 

...oslc fea^e. ^ . " „f .oric, Mza^dous. o. types of »as«. 

i^ventionalsohasapphcafoaforclea^ ' ^ed to visualize a.d ia.e,pr.t data 
- — - le«a.ea.apa.c..si.. Sue. 

'^'^^""^«Cr:llrp.o..i..S0leanupa.v 
visualization and mterpretation ib u 

d.«aciea.-uppianf<.ap^e....e^^^^ ^^^^^^^^ 
„ rr:;r!Lr— e.e.e....e.o.i..«...o.a 
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. gametic resonance imaging (MRI) 
(co^pu^rized a^ai to„>o^aphy) ^ca^er or a magneuc 

procedure. ,h. present invention will be descrilied in 

By way of explanation and example, resen 

. - ^o+a the 3D volume data set. u i& 
.etail below usin. 3D se,snnc ^^^'^ „f 3D seisn>io 

5 however, « the present invenaonrs no. any 

. ■ rti..ulailvuseit.lasavisualizationtoolforinterpreting 
^epresentinventtonrspa^c^^y-^^ 

'^-""""'r IpresentinSDvol^edatase... Suohphysiealfeantxesor 

10 paraMe.ersthatareni.erenUVP--"" ^j^^^.^;,, 3D data 

p^anteters »e typically not apparent fio» eonven ^ 

— — aosertom^^^^^ 

Be^^ of «s tnental re»r..™etion, it is difflcnl. to 
teee-dimensional teatnre. Because o 3^ 

, — ateandsbareanton^t.^-^^^^^^^^ 

cental image of.be terratn wtU no. " ^ " Throu^tbe 

^..dsoranaly^a^— ^ 

use of 3D computer graphics, users seismic data, a team of geologists, 

^o+ocP-f; Bv visualizmg 3D seismic udto, 
SDintageofSDvolt^edatasets By exploration and 

oeophysieists, and engineers can rnterpretftevrsuai 

20 geopny . ,h »s diillms locations and weU paths. 

pioduoUondects,ons.suchasdrJlmg j^^^tion uses the computer 

. Toaccomplish^e— — 2^ . , 

.apbics technique, of ~ ^^^^ . . . 3D o.ect. .or 

„ean. wrapping or mappmg a 2D ^ ^^^^^ coffee oup. 
,3 ex^ple,apbotographofaper»n.an-^^^ 

^''^'^"■•^reo eys.oaviewerthe.hree— onalna^^ 
dimensionalobjectmamanner.ha.eon y ^^.^^^^ display or 

- o.e. e.n though .e ,.ting, and shadmg 

screen. ^^^^^ ^^^^^ 

techniques to convey to the mmd of the ^^^^ ^^^^^ 

Ughtordarlcetc. THe perspective or viewpomt can be 
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all sides of the 3D object. Volume rendering typically includes some type of 
transparency/opacity (opacity = 1 - transparency) control so that certain parts of the 3D 
object are more transparent, thereby allowing a viewer to "see through" outer surfaces of 
an object and view its internal structures. 

The present invention enables fast visualization and analysis of very large 3D 
volume data sets thi-ough the use of a "sampling probe.", also referred to herein as a 
"probe" or "probe object". As explained in more detail below, the sampling probes of the 
present invention have numerous attributes, one of which is that they are typically created 
as a 3D sub-volume of the whole 3D volume data set to be \asualized and analyzed. 

A number of sampling probes can be created, shaped, sized, and moved 
interactively by the user withui the whole 3D volume data set. The intersection of the 
sampling probe with the whole 3D volume data set is texture mapped onto the surfaces 
of the sampling probe, or volmne rendered with varying degrees of transparency within 
the sampling probe. As used herein, the term "interactive" or "interactively" refers to 
1 5 changing or re-drawmg an image m response to user input at a rate sufficiently fast that 

the user perceives an instantaneous or real-time change in the image, without perceptible 
delay or lag. In practice, a frame rate of at least about 10 to 15 frames per second is 
sufficient to achieve interactive imaging as described herein. For example, as the user 
moves the sampling probe, such as by "clicldng and dragging" with a "mouse", the user 
20 perceives the textm^e on tlie surfaces of the sampling probe changing in "real-time" with 

movement of the sampling probe. As the sampling probe changes shape, size, or location, 
there is no user-perceivable delay or lag in imaging the texture, or, with varying degrees 
of transparency, volume-rendered attributes. In this maimer, the user can interactively 
move the sampling probes through the whole 3D voliame, and more easily and effectively 
25 visualize and interpret the features and physical parameters that are present within the 

geographic space represented by the whole 3D volume data set. 

System Description 

The present invention may be implemented using hardware, software or a 
3 0 combination thereof, and may be implemented in a computer system or other processing 

system. One embodiment of a software or program structure 1 GO for implementing the 
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present invention is shown in FIG. 1. At the base of program structure 100 is an 
operating system 1 02. Suitable operating systems 102 include, for example, the UNIX® 
operating system, or Windows NT® from Microsoft Corporation, or other operating 
systems as would be apparent to one of skill in the relevant art. 
5 Menu and vraidowing software 104 overlays operating system 102. Menu and 

windowing software 104 is used to provide various menus and windows to facilitate 
interaction with the user, and to obtain nser input and instructions. Menu and windowing 
software 104 can include, for example, Microsoft Windows™, X Window System™ 
(registered ti*ademark of Massachusetts Institute of Technology), and MOTIF™ 

1 0 (registered trademark of Open Software Foundation Inc.). As would be readily apparent 

to one of skill in the relevant art, other menu aad windowing softvs^are could also be used. 

A basic graphics library 1 06 overlays menu and windowing software 1 04. Basic 
graphics library' 106 is an application programming interface (API) for 3D computer 
graphics. The fimctions performed by basic graphics library 106 include, for example, 

15 geometric and raster primitives, RGB A or color index mode, display list or immediate 

mode, viewuig and modeling transformations, lighting and shading, hidden surface 
removal, alpha blending (translucency), anti-aliasmg, texture mapping, atmospheric 
effects (fog, smoke, haze), feedback and selection, stencil planes, and accumulation 
buffer. 

20 A particularly preferred basic graphics librar^^ 1 06 is OpenGL®, available from 

Silicon Graphics, Inc. ("SGI"), Mountain View, California. The OpenGL® API is a 
multiplatform industrj^ standard tliat is hardware, window, and operating system 
mdependent. OpenGL® is designed to be callable from C, C-H-, FORTRAN, Ada and 
Java programming languages. OpenGL® performs each of the functions Usted above for 

25 basic graphics library 1 06. Some commands in OpenGL® specify geometric objects to 

be drawn, and others control how the objects are handled. All elements of the OpenGL® 
state, even the contents of the texture memory and the frame buffer, can be obtained by 
a client application using OpenGL®. OpenGL® and the client application may operate 
on the same or different macliines because OpenGL® is network transparent. OpenGL® 

30 is described in more detail in tho OpenGL® Programming Guide (ISBN: 0-20 1 -63274-8) 
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and the OpenGL® Reference Manual (ISBN: 0-201-63276-4), the entirety of both of 
which are incorporated herein by reference. 

Visual simulation graphics library 1 08 overlays basic graphics hbrary 1 06. Visual 
sunulation graphics library 108 is an API for creating real-time, multi-processed 3D 
visual simulation graphics applications. Visual simulation graphics library 1 08 provides 
functions that bundle together graphics library state control functions such as lighting, 
materials, texture, and transparency. These functions track state and the creation of 
display lists that can be rendered later. 

A particularly preferred visual simulation graphics library 1 08 is IRIS Performer, 
available from SGI in Mountain View, California. IRIS Performer supports the 
OpenGL® graphics Ubrary discussed above. IRIS Performer includes two main libraries, 
libpf and libpr, and four associated libraries, hbpfdu, libpfdb, libpflii, and libpfutiL 

The basis of IRIS Performer is tlie performance rendering library libpr, a low- 
level library providing high speed rendering functions based on GeoSets and graphics 
state control using GeoStates. GeoSets are collections of drawable geometry that group 
same-type graphics prunitives (e.g., triangles or quads) into one data object. The GeoSet 
contaiQS no geometry itself, only pointers to data arrays and index an*ays. Because all the 
primitives in a GeoSet are of the same type and have tiie same attributes, rendering of 
most databases is performed at maximimi hardware speed. GeoStates provide graphics 
state definitions (e.g., texture or material) for GeoSets. 

Layered above libpr is libpf, a real-tune visual simulation environment providing 
a high-performance multi-process database rendering system that optimizes use of 
mialtiprocessing hardwai'e. The database utility library, libpfdu, provides functions for 
defining both geometric and appearance attributes of 3D objects, shares state and 
materials, and generates triangle strips from independent polygonal input. The database 
library libpfdb uses the facilities of libpfdu, libpf, and libpr to import database files in 
a number of industry standard database formats. The libpfui is a user interface library 
that provides building blocks for writing manipulation components for user interfaces (C 
and C-f-f programming languages). Finally, the Ubpfutil is the utility library that provides 
routines for hnplementing tasks such as Multichannel Option support and graphical user 
interface (GUI) tools. 
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that uses IRIS Performer and OpenGL® API typrcally 
A. application program that .ses IRI ^^^^^ 

carry out the following steps ia preparmg for reai-Ume 3D 
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25 



3. 
4. 
5 
6. 
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Initialize IRIS Performer; 

specify number of graphics pipelines, choose 
^.Itiprocessing configuration, and specify hardware mode as 

needed; 

Initialize chosen multiprocessing mode; 
Initialize frame rate and set frame-extend policy; 
Create, configure, and open windows as required; and 
, Createandconfiguredisplaychannelsasveqmred 
Oncetheapplicat.onprogran.hascreatedagraphicalrendermgenv.orn^^^^^^^ 

nig o" steps 1 through 6 ahove, .en the application program typrcally 
iterates thr-ough a main simulation loop once per fremre. 
7. Compute dynamics, update model matrices, etc.; 
g Delay until the next frame time; 

9. Perform latency critical viewpoint updates; 

10. Draw a frame. 

carried out by, eacn oi vibu i niin a manner Icnown 

. ^ ■ rr c^fHvare 104 and operatmg system 102 m a mami 
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referred to herein as sampling probes or probes. Program 
objects so that they have the following attributes. 
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A probe corresponds to a sub-volume of a larger 3D volume. Particularly, a probe 
defines a sub-set that is less than the complete data set of voxels for a 3D volume data set. 
A probe could be configured to be equal to or coextensive with the complete data set of 
voxels for a 3D voliraie data set, but the functionality of the present invention is best 
carried out when the probe corresponds to a sub-volume and defines a sub-set tliat is less 
than the complete data set of voxels for a 3D volume data set. For example, a 3D volume 
data set of seismic data can contain fi'om about 500 (megabytes) to about 10 GB 
(gigabj^es) or more of data. A 2,500 square kilometer geographic space of typical 3D 
seismic data contains about 8 GB of data. A probe of the present invention for a 500 MB 
seismic data set would preferably contain about 10-20 MB of data. 

By using probes that are a sub-volume of the larger 3D volume, the quantity of 
data that must be processed and re-drawn for each frame of an image is dramatically 
reduced, thereby increasing the speed with which the image can be re-drawn. The 
volxmie of a three-dimensional cube is proportional to the tliird power or "cube" of the 
dimensions of the three-dunensional cube. Likewise, the quantity of data in a 3D volume 
data set is proportional to the third power or "cube" of its size. Therefore, tlie quantity 
of data in a sub-volume of a larger 3D volume will be proportional to the "cubed root" 
(V) of the quantity of data in the larger 3D volume. As such, the quantity of data in a 
probe of the present invention will be proportional to the "cubed root" (V) of tlie quantity 
of data in the 3D volume of which it is a sub-volume. By only having to process the sub- 
set of data that relates to the sub-volume of the probe, the present invention can re-draw 
an image in response to user input at a rate sufficiently fast that tlie user perceives an 
instantaneous or real-time change in the image, without perceptible delay or lag. 

The probes of the present invention can be interactively changed in shape and/or 
size, and interactively moved within the larger 3D volume. The outside geometry or 
surfaces of a probe can be interactively drawn opaque or texture mapped while the probe 
is being changed in shape and/or size or while the probe is being moved. The probe can 
be drawn or volume rendered with varying degrees of transparency while the probe is 
being changed in shape and/or size or moved, thereby revealing the internal structures or 
features of the probe. 
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The 3D sampling probes of the present invention can have any shape, including 
rectangular shapes having one or more right angles and non-rectangular shapes having 



no right angles. The 3D sampling probes of the present invention can have orthogonal or 
perpendicular planes as outer surfaces (e.g., squares and rectangles), parallel planes as 
outer surfaces (e.g., parallelograms), or curved outer surfaces (e.g., spheres, ovals, or 
cylinders). The present invention is not limited to 3D sampling probes of any particular 
shape. The 3D sampling probes of the present invention can have arbitrary shapes, such 
as the shape of a geologic feature identified by a user. For example, as a user moves tlie 
3D sampling probe through a 3D volume of seismic data, a geologic feature may be 
visualized and identified by the user. The 3D sampling probe can be interactively re- 
shaped by tlie user to match the shape of the geologic feature, thereby enabling the user 
to better visualize and define the extent of that geologic feature. 

A probe can be used to cut into another probe, and the intersection of the two 
probes can be unaged. A probe can be used to highlight data in accordance with a seed 
selection algorithm. A probe can also be used to "erase" or delete data in accordance with 
a seed de-selection algorithm. These attiibutes will be explained in more detail below. 

FIG. 2 shows a block diagram of one embodiment of 3D sampling probe program 
110. Program 110 includes a User Interface Module (UIM) 210, a Graphics Processing 
Module (GPM) 220, and a Volume Sampling Module (VSM) 230. A 3D volume data 
set is illustrated as data volume 240, also referred to herein as a 3D volume. UIM 210 
and GPM 220 communicate via a bi-directional pathway 212. GPM 220 sends 
iastructions and requests for data to VSM 230 via pathway 222. UIM 210 sends 
instructions and requests to VSM 230 via pathway 214. UIM 210 interacts with data 
volume 240 through pathway 216. 

Voxel data from data volume 240 is transferred to VSM 230 via data pathway 
234. VSM 230 transfers data to GPM 220 via data pathway 232. Data volume 240 stores 
the 3D volume data set in a manner well Imown to one of skill in the relevant art. For 
example, the format for data volume 240 can consist of two parts, a volume header 
followed by the body of data that is as long as the size of the data set. The volume header 
typically contains information, in a prescribed sequence, such as the file path (location) 
of the data set, size, dimensions in the x, y, and z directions, annotations for the x, y, and 
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z axes, annotations for the datavalue, etc. The body of data is a binary sequence of bytes,, 
one or more bytes per data value, tliat can be ordered in the following manner. The first 
byte is tlae datavalue at volume location (x, y, z)=(0,0,0). The second byte is the 
datavalue at volume location (1,0,0), the third byte is the datavalue at volume location 
5 (2,0,0), etc. When the x dimension is exhausted, then the y dimension is incremented, 

and finally the z dimension is incremented. The present invention is not limited in any 
way to a particular data format for data volume 240. 

User Interface Module 210 handles the user interface to receive commands, 
instructions, and input data from the user. UIM 210 interfaces with the user through a 

1 0 variety of menus through which the user can select various options and settings, either 

through keyboard selection or through one or more user-manipulated input devices, such 
as a "mouse", or a 3D pointing device. UIM 210 receives user input as the user 
manipulates the input device to move, size, shape, etc. a 3D sampling probe. 

The primary functions carried out by UIM 210 will now be described. UIM 210 

1 5 inputs from the user the identification of one or more 3D volume data sets (represented 

by data volume 240) to use for imaging and analysis. When a plurality of data volumes 
are used, the datavalue for each of the plurality of data volvimes represents a different 
physical parameter or attribute for the same geographic space. By way of example, a 
plm-ality of data volumes could mclude a geology volume, a temperature volume, and a 

20 water-saturation volume. The voxels in the geology volume can be expressed in the form 

(x, y, 2, seismic amplitude). Tlie voxels in the temperature volume can be expressed in 
the fomi (x, y, z, °C). 

The voxels in the water-saturation volume can be expressed in the form (x, y, z, 
%saturation). The physical or geographic space defined by the voxels in each of these 

25 volumes is the same. However, for any specific spatial location (xq, y©, Zq), the seismic 

amplitude would be contained in the geology volume, the temperature in the temperature 
volume, and the water-saturation in the water-saturation volume. 

UTM 210 inputs firom the user ioformation to create one or more 3D sampling 
probes. Such information includes size, shape, and initial location of the probe. Such 

30 information can also include imaging attributes such as color, lighting, shaduig, and 

transparency (or opacity). By adjusting opacity as a function of datavalue, certain 
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portions of the data volume are more ti*ansparent, thereby allowing a viewer to see 
through surfaces. An exemplary opacity curve 300 is shown in FIG. 3. Opacity curve 
300 illustrates opacity (1 -transparency) as a function of datavalue. As would be readily 
apparent to one skilled in the art, datavalues with greater opacity (less transparency) will 
mask the imaging or display of datavalues with lower opacity (more transparency). 
Conversely, datavalues will less opacity and greater transparency will permit the imaging 
or display of datavalues with greater opacity and lower transpai'ency. 

UIM 210 receives input from the user for sizing and shaping the 3D sampling 
probes. As described in more detail below, in a preferred embodiment of the present 
invention, the user changes the shape and/or size of aprobe by cliclcing onto "sizing tabs" 
on the probe, and maldng changes in the dimensions of the probe in one or more 
dhections. UIM 210 receives input from the user to move the position or location of a 
3D sampling probe within the data volume. In a preferred embodiment, a user 
manipulates a mouse to "click" onto a surface of the probe to be moved, and then moves 
the mouse to move the probe throughout the geographic space defined by tlie data 
volume. 

UIM 210 receives input from the user to carry out "auto picking" processes. In 
an auto picking process, data points (voxels) are selected based upon a selection 
algorithm. In a preferred embodiment, the selection algorithm is based upon a seed point 
within the 3D data volume. The selection algorithm then selects data points that: (i) 
satisfy the selection criteria or algorithm (e.g.^ have a datavalue within a specified filter 
range); and (ii) have a connectivity witli or are comiected to the seed point. Through 
UIM 210, the user is prompted to identify a seed point within the 3D volvime, and to 
identify a filter range of datavalues used by the selection algorithm to "pick" the selected 
points. Preferably, the seed point is within one of the 3D sampling probes. 

UIM 210 also receives input from tlie user regarding the content of the displayed 
image. For example, the user can preferably select the content of the displayed image. 
The content of the displayed image coxild include only the 3D sampling probe, i.e., its 
intersection with the 3D volume. Additionally, the 3D sampling probe could be 
displayed either with or without a bounding box that defines the outer geometry of the 
probe. Alternatively, the displayed image could include the 3D sampling probe, as well 
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as the data that occupies the background xz, yz, and x}^ planes, and/or the data that 
occupies the 3D volume outside of the 3D sampling probe(s) being displayed. 

To carry out the foregoing functions, UIM 210 sends a request to Voliime 
Sampling Module 230 to load or attach those 3D volume data sets identified by the user. 
UIM 210 communicates via pathway 212 with Graphics Processing Modiale 220 that 
carries out the display and. imaging. 

The primary functions carried out by GPM 220 will now be described. GPM 220 
processes data for imaging of 3D sampling probes with the hghting, shading, 

transparency, and other attributes selected by the user. To do so, GPM 220 uses the 
fiznctions available through basic grapliics library 106 and visual simulation graphics 
library 108 described above. The user can select (tlirough UIM 210) to display only the 
one or more 3D sampling probes that have been created. Alternatively, the user can 
select to display one or more 3D samplmg probes, as well as the 3D data volume outside 
of the probes, i.e. voxels v^thin the 3D volxmie that do not intersect any of the 3D 
sampling probes that ai-e being displayed. 3D sampling probes that are being displayed 
are referred to herein as "active probes". 

GPM 220 processes the re-shaping and move requests that are received by UIM 
210 fi-om the user. GPM 220 draws the re-shaped 3D sampling probe in accordance with 
the user-selected attributes (color, lighting, shading, transparency, etc.). As the user 
inputs a change in shape for a 3D sampling probe, the image with selected attributes is 
re-drawn sufficiently fast to be perceived as real-time by the user, Similai'ly, GPM 220 
draws the 3D sampling probe in the new position or location in accordance with the user- 
selected attributes (color, lighting, shading, transparency, etc.). As the user moves the 
3D sampling probe tlirough the 3D volume, the image of the 3D sampling probe with 
selected attributes is re-drawn sufficiently fast to be perceived as real-time by the user. 

GPM 220 processes "auto piclcing" requests that are received by UIM 210. GPM 
220 will image selected points withui the 3D volume in accordance with the selection 
algorithm. Altematively, GPM 220 will "erase" selected points within the 3D volume 
in accordance with the selection algorithm. 

To carry out the foregomg functions, GPM 220 commxmicates via pathway 212 
with UIM 210 so that the information requested by the user is imaged or displayed with 
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up Step 402 can be carried out so that the data that occupies the background xz, yz, and 
xy planes, or the data that occupies the 3D volume outside of the default 3D sampling 
probe, is not displayed or imaged. Preferably, the present invention is carried out so that 
the user can selectively display, or not display, the data that occupies the background xz, 
5 yz, and xy planes, as vv^ell as the data that occupies the 3D volume outside of the active 

probes. 

In a step 404, UIM 210 is waiting to respond to user input or request. User input 
is received through a user input device suitable for use with a computer, including but not 
limited to, a keyboard, mouse, joystick, ti'ackball, roUerbali, roller point, or other type of 

10 suitable pointing device, etc. Preferably, the user input device comprises a mouse or 

other similar device that enables the user to "click" on a particular displayed image, and 
"drag" that displayed image to another location. Such a user input device allows a user 
to move and re-shape displayed probes. Such a user input device also allows, a user to 
activate drop-down menus, and to select the various options for the color, shading, 

15 lighting, and transparency attributes. A keyboard can also be used for entering 

information relating to the selected attributes. 

Reference numeral 406 refers generally to a plurality of functions that can be 
carried out by the present invention. These functions can be carried out individually or 
simultaneously, depending upon input from the user. For example, a probe can be moved 

20 (function 430) and rotated (function 450) simultaneously. While the functions identified 

by reference numeral 406 are being carried out, the image of the 3D sampling probes is 
being re-drawn sufficiently fast to be perceived as real-time by the user. Each of the 
functions identified by reference numeral 406 will now be described. 

If a user wants to change the default probe, then function 410 is carried out. The 

25 steps for carrying out function 410 are shown in FIG. 5 by way of flow diagram 

corrector 5 A. In a step 502, the changes to the default probe are input by UIM 210 from 
the user. For example, the changes to the default probe can be to the shape or size, the 
location, or the attributes such as color, shading, lighting, and transparency. 

In a step 504, UIM 210 sends a request to GPM 220 to draw the changed default 

30 probe. In a step 506, GPM 220 requests data for the changed default probe from VSM 

230. In making this request, GPM 220 would invoke function 430 if it was necessary to 
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move the default probe, function 440 to re-shape the default probe, and functions 450 or 
460 to rotate the default probe. The foregoing functions will be described in more detail 
below. 

The data that will be extracted from data volume 240 by VSM 230 in response 
5 to the request made by GPM 220 in step 506 will depend upon attributes that have been 

selected by the user. If the opacity settings selected by the user are such that all 
datavalues are opaque, then the data extracted by VSM 230 will be limited to the surfaces 
of the changed default probe. Because of the selected opacity, it will not be possible for 
the user to see inside the changed default probe, so only the data corresponding to the 

1 0 surfaces or outside of the changed default probe will be extracted by VSM 230. In a step 

508, GPM 220 processes the data extracted by VSM 230 for the surfaces of the changed 
defa\ilt probe, and draws the changed default probe by texture mapping onto the surfaces 
in accordance with the attributes selected by the user. By extracting only the data that 
can be seen by the user, the image of the changed default probe can be drawn more 

15 quicldy because less data needs to be processed, i.e., the data corresponding to the 

"inside" of the changed default probe is not processed. 

Alternatively, if the opacity settings selected by the user are such that some of the 
datavalues are opaque and some of the datavalues are transparent, then the data extracted 
by VSM 230 will include the data corresponding to the entire volvmie of the changed 

20 default probe. Because of the selected opacity and transparency, it will be possible for the 

user to see inside the changed default probe, so data corresponding to the entire volume 
of the changed default probe vAll be extracted by VSM 230. In such a situation, GPM 
220 processes the data extracted by VSM 230 m step 508, and draws the changed default 
probe by volume rendering in accordance with the attiibutes selected by the user. 

25 If a user wants to create additional probes, then function 420 is carried out. The 

present invention is not limited to any particular number of active probes. The steps for 
carrying out function 420 are shown in FIG. 6 by way of flow diagram connector 6A. In 
a step 602, the shape, size, location, attributes, etc. for the additional probes are input by 
UIM 210 from the user. In a step 604, UIM 210 sends a request to GPM 220 to draw the 

30 additional probes. 
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In a step 606, GPM 220 requests data for the additional probes from VSM 230. 
In a manner similar to that described above for changing the default probe, the data that 
is extracted from 3D or data volume 240 by VSM 230 will depend upon tlie opacity 
selected by the user for the additional probes. If the opacity settings selected by the ixser 
are such that all datavalues for the additional probes are opaque, then the data extracted 
by VSM 230 will be limited to the surfaces of the additional probes. Alternatively, if the 
opacity settmgs selected by tlie user for the additional probes are such that some of the 
datavalues are opaque and some of the datavalues are transparent, then the data extracted 
by VSM 230 will include the data corresponding to the entire volumes of the additional 
probes. In this manner, the additional probes can be drawn more quickly by minimizing 
the quantity of data that must be processed. 

In a step 608, GPM 220 processes the data extracted by VSM 230 for the 
additional probes, and draws the additional probes in accordance with the attributes 
selected by the user, either by textm*e mapping onto the surfaces of the additional probes, 
or by volume rendering the entire volimie of the additional probes. 

If a user wants to move a probe, then function 430 is carried out. The steps for 
carrying out function 430 are shown in FIG. 7 by way of flow diagram connector 7A. In 
a step 702, tlie new location for the probe is iaput by UIM 210 from the user. In a 
preferred embodiment of the present invention, the user inputs the new location of the 
probe by clicking a mouse or other type of suitable user input device to snap a pointer 
onto a surface of the probe to be moved. The user changes the location of the probe by 
moving the mouse or other suitable user input device in any direction, thereby dragging 
the probe along a trajectory. 

In a step 704, UIM 210 sends a move request to GPM 220 to draw the probe at 
the new location. GPM 220 requests data for the new location of the probe from VSM 
230. In a manner similar to that described above, the data that is extracted from data 
volume 240 by VSM 230 will depend upon the opacity selected by the user for the probe 
being moved. If the opacity settings selected by the user are such that all datavalues for 
the probe being moved are opaque, then the data extracted by VSM 230 will be limited 
to the siorfaces of the probe being moved. Alternatively, if the opacity settings selected 
by the user for the probe being moved are such that some of the datavalues are opaque 
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and some of tlie datavalues ai*e transparent, then the data extracted by VSM 230 will 
include the data corresponding to the entire volume of the probe being moved. In this 
manner, the probe can be drawn at its new location more quickly by minimizing the 
quantity of data that must be processed. 
5 In a step 708, GPM 220 processes the data extracted by VSM 230 for the probe 

being moved, and dmws the probe at its new location in accordance with the attributes 
selected by the user, either by texture mapping onto the smfaces of the probe being 
moved, or by volxmie renderuig the entire volume of the probe being moved. 

As the user moves the probe, for each new location of the probe, steps 702 

1 0 tlrough 708 are repeated at a rate sufficiently fast that the user perceives tlie image of the 

probe, with texture mapping or volume rendering as appropriate, changing in "real-time*' 
with movement of the probe. The image is being re-drawn at a frame rate sufficiently 
fast to be perceived as real-time by the user. 

If a user wants to re-shape a probe, then function 440 is carried out. As used 

1 5 herein, the term "re-shape" refers to any change in dimension of a 3D sampling probe in 

any dnection. The shape of a 3D sampling probe can be changed, or re-shaped, for 
example, by changing the size in one or more directions, such as by changing a square 
probe into a rectangular probe by increasing the size of the probe in the x direction, and 
decreasing the size of the probe in the y direction. As another example, the shape of a 

20 3D sampling probe can be changed by changing the shape from spherical to rectangular. 

As yet another example, a square 3D sampling probe (equal dimensions in the x, y, and 
z directions) can be re-shaped in accordance with the present invention to be a larger or 
smaller square-shaped probe by changing the size equally in each of the x, y, and z 
directions. The re-shaped probe also has a square shape, but as a larger or smaller square. 

25 The steps for carrying out function 440 are shown in FIG. 8 by way of flow 

diagram connector 8A. In a step 802, the new shape and/or size for the probe is input by 
UIM 210 from the user. In a preferred embodiment of the present invention, the user 
inputs the new shape of a probe by clicking a mouse or other type of suitable user input 
device to snap apomter onto a "sizing tab" of the probe to be re-shaped. As used herein, 

30 a "sizing tab" refers to a designated area on a smface of the probe. Such a designated 

area is preferably displayed in a color that is dijfferent from the colors being used to 
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display the features or physical parameters of the 3D volume data set. When the pointer 
is snapped to the sizing tab, manipulation of the mouse or user input device changes the 
dimensions or proportions of the surface on which the sizing tab is located. When the 
desired size or shape is reached, the user again clicks the mouse or user input device to 
5 release the pointer from the sizing tab. Sizing tabs are illustrated in FIGS. 1 5 and 1 6. The 

sizing tabs are the small dark squares that appear on the surfaces of the probes, along the 
bounding geometry of the probes. The location of the sizing tabs is not limited to the 
bounding geometrj^ of the probes. The user changes the shape of the probe by clickuig 
the mouse or other suitable user input device onto a sizing tab, moving the mouse until 

1 0 the surface being changed has the desired shape, and then releasing the mouse from tlie 

sizing tab. This process can be repeated, if necessar}^ using other sizing tabs on the 
probe until the probe is re-shaped to the desired shape. 

It would be readily apparent to one of sidll in the relevant art how to implement 
such a sizing tab for re-shaping the probes of the present invention. It is to be 

1 5 understood, however, that the present invention is not limited to the use of sizing tabs for 

re-shaping probes, and other suitable methods can be used. For example, the user could 
select from a number of pre-set shapes (e.g., squares, rectangle, cylinders, spheres) by 
activating a drop-down menu, or by scrolling through the shapes by repeatedly clicldng 
a mouse, 

20 In a step 804, UIM 210 sends a re-shape request to GPM 220 to draw the re- 

shaped probe. In a step 806, it is determined whether more data is needed to draw the re- 
shaped probe. For example, if the re-shaped probe is of a shape and size that "fits inside" 
the existing probe, then no more data is needed, and processing continues at a step 810. 
Alternatively, if the re-shaped probe is of a shape and size that falls at least partially 

25 outside of the existing probe, then, in a step 808, GPM 220 requests the data needed for 

the re-shaped probe from VSM 230. In a manner similai'to that described above, the data 
that is extracted from 3D or data volume 240 by VSM 230 will depend upon the opacity 
selected by the user for the probe being re-shaped. If the opacity settings selected by the 
user are such that all datavalues for the probe being re-shaped are opaque, then the data 

30 extracted by VSM 230 will be limited to the surfaces of the probe being re-shaped. 

Alternatively, if the opacity settings selected by the user for the probe being re-shaped 
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are such that some of tlie datavalues are opaque and some of the datavalues are 
transparent, then the data extracted by VSM 230 will include tlie data corresponding to 
the entire volume of the probe beuig re-shaped. In tliis manner, the probe can be drawn 
wdth its new shape more quiclcly by minimizing the quantity of data that must be 
processed. 

In step 810, GPM 220 processes the data extracted by VSM 230 for the probe 
being re-shaped, and draws the probe with its new shape in accordance witli the attributes 
selected by the user, either by texture mapping onto the surfaces of the probe being re- 
shaped, or by volume rendering the entire volume of the probe being re-shaped. 

As the user changes the shape of the probe, steps 802 through 810 are repeated 
at a rate sufficiently fast that the user perceives the image of the probe, with texture 
mapping or volmne rendering as appropriate, changing in *'real-time" with the changing 
sliape of the probe. The image is being re-drawn at a frame rate sufficiently fast to be 
perceived as real-time by the user. 

If a user wants to rotate a probe in 3D space, then function 450 is carried out. In 
fimction 450, the 3D orientation, which is the same for both the 3D volume and the 
probe, is changed, thereby rotating the 3D volume and the probe in space. The steps for 
carrying out fimction 450 are shown in FIG. 9 by way of flow diagram connector 9A. In 
a step 902, the new 3D orientation for the 3D volume and the probe is input by UIM 210 
from the user. In a preferred embodiment of the present invention, the user inputs the 
new orientation by clicking a mouse or other type of suitable user input device to snap 
a pointer onto an axis of the probe to be rotated. Manipulation of tlie mouse or user input 
device changes the orientation of that axis. When the desired orientation is reached, the 
user again clicks the mouse or user input device to release the pointer from the axis. It 
would be readily apparent to one of skill in the relevant art how to implement such a 
change in orientation. It is to be understood, however, that the present invention is not 
limited to changing the orientation in this manner. For example, the user could select 
from a number of pre-set rotations (e.g., rotate 90' to the left or right; rotate 45" to the left 
or right, etc.) by activating a drop-down menu, or by scrolling through the rotations by 
repeatedly clicking a mouse. 
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In a step 904, UIM 210 sends a request to rotate in 3D space to GPM 220 to draw 
the rotated probe. In a step 906, GPM 220 requests data for the rotated probe ifrom VSM 
230. In a manner siroiiar to that described above, the data that is extracted from 3D or 
data volume 240 by VSM 230 will depend upon the opacity selected by the user for the 
5 probe being rotated. If the opacity settings selected by the user are such that all 

datavalues for the probe being rotated are opaque, then the data extracted by VSM 230 
will be limited to the surfaces of the probe being rotated. Alternatively, if the opacity 
settings selected by the user for the probe being rotated are such that some of the 
datavalues are opaque and some of the datavalues are transparent, then the data extracted 

1 0 by VSM 230 will include the data correspondmg to the entire volume of the probe being 

rotated. In this manner, the probe can be drawn with its new orientation more quickly by 
minimizing the quantity of data that must be processed. 

In step 908, GPM 220 processes the data extracted by VSM 230 for the probe 
being rotated, and draws the probe with its new orientation in accordance with the 

1 5 attributes selected by the user, eitlier by texture mapping onto the surfaces of the probe 

being rotated, or by volume rendering the entire volume of the probe being rotated. 

As the user rotates the probe in 3D space, steps 902 through 908 are repeated at 
a rate sufficiently fast that the user perceives the image of the probe, with texture 
mapping or volume rendering as appropriate, changing m "real-time" with the changing 

20 orientation of the probe. The image is being re-drawn at a frame rate sufficiently fast to 

be perceived as real-time by the user. 

If a user wants to rotate a probe while it is fixed in 3D space, then function 460 
is carried out. In function 460, the 3D orientation of the probe is rotated independently 
of the 3D orientation of the 3D volume, thereby rotatbig the probe while it is fixed in the 

25 3D space defined by the orientation of the 3D volume. In this manner, the background 

planes for an active probe can be displayed in a fixed orientation, and the active probe 
can be rotated within the background planes. 

The steps for carrying out function 460 are shown in FIG. 10 by way of flow 
diagram connector 1 OA. In a step 1002, the new 3D orientation for the probe is input by 

30 UIM 210 from the user. In a preferred embodiment of the present invention, the user 

selects the option to rotate while fixed in space, for example, from a "drop-down" menu. 



WO 02/37422 




-28- 



PCT/USOO/29835 

The user then hiputs the new orientation for the probe by clickmg a mouse or other type 
of suitable user input device to snap a pointer onto an axis of the probe to be rotated. 
Maxiipulatiou of the mouse or user input device changes the orientation of that axis. 
When the desired orientation is reached, the user again chcks the mouse or user input 
5 device to release the pointer from the axis. It would be readily apparent to one of skill 

in the relevant art how to implement such a change in orientation. It is to be understood, 
however, that the present invention is not limited to changing the orientation in this 
mamier. For example, the user could select from a number of pre-set rotations (e.g., 
rotate 90*" to the left or right; rotate 45' to the left or right, etc.) by activating a drop-down 

10 menu, or by scrolling through the rotations by repeatedly clicking a mouse. 

In a step 1004^ UIM 210 sends a request to rotate wliile fixed in space to GPM 
220 to draw tire rotated probe. In a step 1006, GPM 220 requests data for the rotated 
probe from VSM 230. In a maimer similar to that described above, tlie data that is 
extracted from 3D or data volume 240 by VSM 230 will depend upon the opacity 

15 selected by the user for the probe being rotated. If the opacity settings selected by the 

user are such that all datavalues for tlie probe being rotated are opaque, then the data 
extracted by VSM 230 will be limited to the surfaces of the probe being rotated. 
Alternatively, if the opacity settmgs selected by the user for the probe being rotated are 
such that some of the datavalues are opaque and some of the datavalues are transparent, 

20 then the data extracted by VSM 230 will include the data corresponding to the entire 

volume of the probe being rotated. In tliis manner, the probe can be drawn with its new 
orientation more quickly by minimizing the quantity of data that must be processed. 

In step 1008, GPM 220 processes the data extracted by VSM 230 for the probe 
being rotated, and draws the probe with its new orientation m accordance with the 

25 attributes selected by the user, either by texture mapping onto the surfaces of tlie probe 

being rotated, or by volume rendering the entire volume of the probe being rotated. 

As the user rotates the probe while it is fixed in space, steps 1002 through 1008 
are repeated at a rate sufficiently fast that the user perceives the image of the probe, with 
texture mapping or volume rendering as appropriate, changing in "real-time" with the 

30 changing orientation of the probe. The image is being re-drawn at a frame rate 

sufficiently fast to be perceived as real-time by the user. 
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If a user wants to carr)^ out an "auto picking" process, then function 470 is carried 
out. The steps for carrying out function 470 are shown in FIG. 11 by way of flow diagram 
connector 1 1 A. In a step 1102, a seed pomt within the data set of the 3D volume, and a 
selection criteria based on datavalues, are input by UIM 210 from the user. Preferably 
tlie seed point is within the data set of voxels that defines a probe. As described below, 
such a probe is referred to herein as a seed 3D sampling probe or an eraser 3D sampling 
probe. However, the seed point can be witliin the data set of voxels that defines the 3D 
volume, outside of an active probe. In a preferred embodiment of the present invention, 
tlie user selects the option to execute an auto picking process, for example, from a "di'op- 
down" menu. The user then selects the seed point by clicking a mouse or other type of 
suitable user input device to snap a pointer onto the desired seed point. The selection 
criteria can be input, for example, by graphically selecting a range, or by keying in 
specific nimierical values. It would be readily apparent to one of sldll in the relevant art 
how to input from the user a seed point and filter range of datavalues. 

In a step 1104, UIM 210 sends an auto piclcing request to GPM 220 to draw the 
rotated probe. In a step 1106, GPM 220 requests selected points to image from VSM 230. 
The selected points are those that are connected to the seed point, and that have a 
datavalue within the selection criteria. 

In step 1108, GPM 220 processes the data extracted by VSM 230 to draw the 
selected points. The selected points are preferably highlighted by being drawn in a color 
different from those used to depict the features or physical parameters of the 3D volume 
data set. Alternatively, step 1108 can be carried out to "erase" or delete from the image 
the selected points. 

In a similar manner, auto picking fiinction 470 can be used to "erase" or de-select 
points. For example, an eraser 3D sampling probe is defined, such as by invoicing 
fiinction 420 to create an additional probe. A "de-selection" criteria based on datavalues 
is defined. Points previously selected by an auto picking operation that satisfy the de- 
selection criteria are identified as candidates for de-selection. As the eraser 3D sampling 
probe moves through the 3D volume, the de-selected points are deleted from the image, 
and the image is re-drawn sufficiently fast to be perceived as real-time by the user. 
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Once auto picking function 470 is initiated by the user, it can be carried out 
simultaneous!)^ with, for example, move function 430. In this manner, as the user moves 
the probe, steps 1102 through 1108 (and steps 702 through 708) are repeated at a rate 
sufficiently fast that the user perceives the image of the probe, with the selected points, 
changmg m "real-time" with the changing location of the probe. As the probe is moved, 
the selected points can be highlighted by being drawn in a suitable color, thereby having 
the auto-picldng 3D sampling probe function as a "highlighter" as it moves through the 
3D volume. Alternatively, as the probe is moved, points previously selected by an auto 
picldng operation can be "erased'* or deleted &om the image, thereby having the probe 
function as an "eraser" or eraser 3D sampling probe as it moves through the 3D volume. 
In either embodiment, the image is being re-drawn at a frame rate sufficiently fast to be 
perceived as real-time by the user. 

If a user wants to create a "ribbon section," then function 480 is performed. The 
steps necessary for performing function 480 are described further below in reference to 
FIG. 17 and by way of block diagram connector 18A in FIG. 18. 

If a user wants to create a "3D surface' representative of a physical phenomena 
found witliin a 3D volume data set, then function 490 is performed. The steps necessary 
for performing function 490 are described further below in reference to FIG. 19 and by 
way of block diagram connector 20 A in FIG. 20. 

In any event where a user desires to carry out one or more of the functions 
described above such as more probe (430), re-shape probe (440), create a ribbon section 
(480) and create a 3D surface (490), each function can be performed independent of, or 
in connection with, one or more of the other functions. 

With reference now to FIG. 12, one embodiment of a computer system suitable 
for use with the present invention is shovm. A graphics supercomputer 1210 contains one 
or more central processing imits (CPU) or processors 1212. Supercomputer 1210 
contains arandom access memory (RA]\4) 1214 that can be accessed by processors 1212. 
Supercomputer 1210 also contains one or more graphics modules 1216 that also access 
RAM 1214. Graphics modules 1216 execute the functions carried out by Graphics 
Processing Module 220, using hardware (such as specialized graphics processors) or a 
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combination of hardware and sofhvai*e. A user input device 1218 allows anser to control 
and input information to grapliics supercomputer 1210. 

A particularly preferred graphics supercomputer is an Onyx2 Infinite Reality 
system, available &om Silicon Graphics, Inc., Mountain View, CA, configured with eight 
5 processors, three graphics pipelines, 16 GB of main memory, and 250 GB of disk 

memory. Such a graphics supercomputer has a scalable, high-bandwidth, low-latency 
architecture to provide high speed rendering on multiple graphics pipelines. Graphics 
supercomputers from other vendors, such as Hewlett-Packard Company of Palo Alto, CA 
or Sun Microsystems of Mountain View, CA could also be used. 

10 The graphics data forming the image to be displayed is sent fi*om graphics 

supercomputer 1210 to a multiple-screen display system 1220 for projection onto a 
screen 1230. In the embodiment shown in FIG. 12, three projectors are used. From the 
perspective of a user viewing the image on screen 1230, the three proj ectors include a left 
projector 1240, a center projector 1250, and a right projector 1260. Although three 

15 projectors are shown, the present invention is not liirdted to the use of any particular 

number of projectors. 

Projector 1240 has a projection field on screen 1230, shown generally at 1242, 
between a point 1241 and a point 1243. Projector 1250 has a projection field on screen 
1230, shovvoi generally at 1252, between a point 1251 and a point 1253. Projector 1260 

20 has a projection field on screen 1230, shown generally at 1262, between a point 1261 and 

a point 1263. Projection fields 1242 and 1252 have an overlap region 1244, between 
points 1251 and 1243. Similarly, projection fields 1262 and 1252 have an overlap region 
1264, between points 1261 and 1253. The image to be displayed is divided into three 
(left, center, and right) over-lapping sub-images. By simultaneously projecting the three 

25 over-lapping sub-images, the field-of-view to the user is increased over that available, for 

example, on a monitor or through the use of just one projector. As an example, use of the 
three over-lapping sub-images shovm in FIG. 12 increases the field-of-view to 
approximately 160°. Overlap regions 1244 and 1264 are each approxhnately 53\ 
Multiple-screen display system 1220 accounts for overlap regions 1244 and 1264 in a 

30 well-lcnown manner to edge-blend the images of the three projectors to form one 
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seamless image on screen 1230. Suitable display and projector systems are available 
from SEOS, London, England, such as the Barco projector units. 

FIG. 1 3 shows an alternate embodiment of a computer system suitable for use 
with the present invention. In the embodiment shown in FIG. 13, graphics 
supercomputer 1210 is configured with multiple processors 1212, RAM 1214, and two 
graphics modules 1216. Graphics workstations suitable foruse in the embodiment shown 
in FIG. 1 3 are available from Silicon Graphics, Inc. or Sun Microsystems, Each graphics 
module 1216 is connected to amonitor 1320 for display. Monitor 1320 should preferably 
be a color graphics monitor suitable for display of graphics such as that shown in FIGS. 
15 and 16. Preferably, one of monitors 1320 displays the image of the 3D sampling 
probes, and the other of monitors 1320 displays the various menus used to operate 3D 
sampluig probe program 1 10. FIG. 13 also shows a keyboard 1330 and a mouse 1332 
that function as user input devices. 

A computer system capable of carrying out the functionality described herein is 
shown in more detail in FIG. 14. Computer system 1402 includes one or more 
processors, such as processor 1404. Processor 1404 is connected to a communication bus 
1406. Various software embodiments are described in terms of this exemplary computer 
system. After reading this description, it will become apparent to a person skilled in the 
relevant art how to implement die invention using other computer systems and/or 
computer architectures. 

Computer system 1402 also includes a main memory 1408, preferably random access 
memory (RAM), and can also include a secondarj^ memory 1410. Secondary memory 
1410 can include, for example, a hard disk drive 1412 and/or a removable storage drive 
1414, representing a floppy disk drive, a magnetic tape drive, an optical disk drive, etc. 
Removable storage drive 1414 reads from and/or writes to a removable storage unit 1418 
in a well known maimer. Removable storage unit 1418, represents a floppy disk, 
magnetic tape, optical disk, etc. which is read by and written to by removable storage 
drive 1414. As will be appreciated, removable storage unit 1418 includes a computer 
usable storage medium having stored therein computer software and/or data. 

In alternative embodiments, secondary memory 1410 may include other similar 
means for allowing computer programs or other instructions to be loaded into computer 
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system 1402. Such means can include, for example, a removable storage unit 1422 and 
an interface 1420. Examples of such can include a program cartridge and cartridge 
interface (such as that found in video game devices), a removable memory chip (such as 
an EPROM, or PROM) aad associated socket, and other removable storage units 1422 
5 and interfaces 1420 which allow software and data to be transferred from removable 

storage unit 1422 to computer system 1402. 

Computer system 1402 can also include a communications interface 1424. 
Communications interface 1424 allows software and data to be transferred between 
computer system 1402 and external devices. Examples of communications interface 

10 1424 can include a modem, a network interface (such as aia Ethernet card), a 

communications port, a PCMCIA slot and card, etc. Software and data transferred via 
communications interface 1424 are in the form of signals 1426 that can be electronic, 
electromagnetic, optical or other signals capable of being received by communications 
interface 1424. Signals 1426 are provided to communications interface via a channel 

15 1428. Channel 1428 carries signals 1426 and can be implemented using wire or cable, 

fiber optics, a phone line, a cellular phone linic, an RF linlc and otlier communications 
channels. 

In this document, the terms "computer program medium" and "computer usable 
medium" are used to generally refer to media such as removable storage device 1418, a 
20 hard disk installed in hard disk drive 1412, and signals 1426. These computer program 

products are meaas for providing software to computer system 1402. 

Computer programs (also called computer control logic) are stored in main 
memory 1 408 and/or secondary memory 1410. Computer programs can also be received 
via communications interface 1424. Such computer programs, when executed, enable 
25 computer system 1402 to perform the features of the present invention as discussed 

herein. In particular, the computer programs, when executed, enable processor 1404 to 
perform the features of the present invention. Accordingly, such computer programs 
represent controllers of computer system 1402. 

In an embodiment where the invention is implemented using software, the 
3 0 software may be stored in a computer program product and loaded into computer system 

1402 using removable storage drive 1414, hard drive 1412 or communications interface 
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1424. The control logic (software), when executed by processor 1404, causes processor 
1404 to perform the functions of the invention as described herein. 

In another embodiment, the invention is implemented primarily in hardware 
• using, for example, hardware components such as application specific integrated circuits 
5 (ASICs). Implementation of such a hardware state machine so as to perform the 

functions described herein will be apparent to persons skilled in the relevant art(s). 

In yet another embodiment, the invention is implemented using a combination of 
both hardware and software. 
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1 0 System Operation and Results 

The operation and results of the present invention will now be described, using 
a data volume 240 that contains seismic data (datavalues representing seismic 
amplitudes). The user specifies the particular seismic data volume to be used, which is 
loaded from disk into main memory. A default 3D sampling probe is drawn. The user 

1 5 specifies the colors to be used for the seismic amphtudes. The degree of transparency can 

also be selected. The three probes shown in FIG. 15 are all opaque, with the intersection 
of the probes and the seismic data volume texture mapped onto the surfaces of the probes. 
One of the probes is displayed with the boimding geometry shown; the other two probes 
axe displayed without the boxmding geometry. 

20 FIG. 1 5 shows three active probes. The user has selected not to display the data 

contained in the backgroimd planes and in the remainder of the seismic data volume 
outside of the active probes. Two of the probes shown in FIG. 15 intersect each other, 
and the intersection of the two probes is displayed. In this manner, the user can more 
readily visualize and interpret geologic features inherent in the seismic data volume. For 

25 example, a geologic feature, represented by a dark band between two light bands, extends 

across the face of the larger intersecting probe and "turns the comer" to extend onto the 
face of the smaller intersecting probe perpendicular to it. The ability to move and 
intersect the probes with each other, and throughout the seismic data volume, enables a 
user to better interpret and track the extent of such a geologic feature. 

30 FIG. 1 6 illustrates how one probe can be used to "cut" another probe to create a 

"hole" in a probe. As with FIG. 1 5, the user has selected not to display the data contained 
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in the background planes and in the lemainder of the seismic data volume outside of the 
active probes. FIG. 1 6 illustrates that the opacity settings can be individually selected by 
the user for each active probe. One of the probes shown in FIG. 16 is opaque, so that it 
is not possible to see through the surfaces of this probe. In order to see internal to this 
outer probe, it must be cut away by another probe. The outer opaque probe is referred 
to as a "data probe". A second completely transpai*ent "cut probe" has been used to cut 
out a 3D sub-section of the data probe. Because the cut probe is completely transparent, 
it is not visible in FIG. 16. However, the fact that the completely transparent cut probe 
is present is evidenced by the fact that the opaque internal surfaces of the data probe are 
visible. The image of the intersection of the data probe and the cut probe is the 
intersecting surface internal to the data probe. 

A third active probe is shown in FIG. 1 6. The tliird probe is displayed with the 
bounding geometry shown. The third probe is volume rendered with varying degrees of 
transparency so that the user can see through the outer surfaces of the probe and view 
geologic features within the third probe. As shown in FIG. 1 6, the third probe is volume 
rendered partially within the 3D sub-section of the data probe that has been cut away by 
the cut probe. 

The thii'd volume-rendered probe shown m FIG. 1 6 also contains selected points 
that have been selected through a seed picking process (function 470). The selected 
points have been imaged in a manner to highlight tiiem for the user . The selected points 
are shown in FIG. 1 6 as connecting points. The seed point is illustrated in FIG. 1 6 by the 
darker sphere. 

FIG. 17 illustrates one embodiment of a ribbon section 1710 in accord with the 
present invention. Like the probe, ribbon sections are a 3 -D volume visualization method 
for displaying data along user-defined traverses through a 3-D volume data set within a 
probe. The traverses cut through the 3-D volume data set somewhat like a cookie cutter 
through dough and therefore are referred to herein as cookie planes such as cookie planes 
1712 and 1714. Ribbon sections may display the 3-D data in orientations that do not 
necessarily conform to the orientation of the ordinate axis of the 3-D volume data set 
and/or probe. Ribbon section 1710 is produced by the user in accord with the process 
described herein below. The traverses or cookie planes are deiSned by digitizing control 
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points such as coatrol points 1716, 1718, and 1720, which may be selected from a face 
of tlie probe referred to as the "probe face plane." The user may create a transparent cut 
probe as described above to interface with an opaque ribbon section. Alternatively, the 
probe can be made opaque and the ribbon section can be made transparent if desired. 
5 The control points may be used to produce a plurality of line segments such as 

line segments 1722, 1724, and 1726, which collectively are referred to as polyline 1728 
which is like a polygon but may or may not be closed. Therefore, tlie Ime sogments may 
form an open or closed line so that a single or multiple cookie planes may be produced. 
In a preferred embodiment, the area of data display, i.e. the cookie planes, is projected 

1 0 along a direction perpendicular to the probe face plane and the data displayed may extend 

to an opposing face of the probe referred to as the "opposing probe face plane." After 
creating a ribbon section, the user may use the mouse controls or keyboard to select, 
move, di-ag, or grab the control points to edit the ribbon section 1710 in real time and 
display different data from the 3D volume data set along tlie cookie planes 1712 and 

15 1714. In addition to editing the ribbon section 1710 to view different data from the 3D 

volume data set vwthin the probe, the entire ribbon section 1710 and probe may be 
simultaneously moved to a different position in order to view different data fi-om the 3D 
volume data set outside the boimdaries of the probe at its prior position. Active control 
point 1720 is preferably high-lighted or colored differently as compared to the other 

20 control points to indicate control point 1720 is in an active state for operations such as 

moving, deleting, or otherwise editing as discussed further below. Control points may 
be inserted or deleted after the initial ribbon section construction. The ribbon section 
geometry and orientation may also be saved for future work sessions. 

FIG. 18 depicts a block diagram of system 1810 for program modules in a 

25 presently preferred embodiment of the invention for producing ribbon sections at real 

time frame rates as discussed herein above. Probe module 1822 provides initial context 
for sketching polyline 1728. Probe module 1822 supplies user activity data such as 
mouse clicks or keystrokes to cookie manager 1824. Thus, user activity data such as 
creation of control points, deletion of control points, moving of control points, moving 

30 of the entire probe, and the like are inserted into system 1810. 
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Coolcie manager 1824 manages the user input data supplied by probe 1822. 
Cookie manager 1824 distributes the data, e.g. control points add, move, and delete as 
appropriate to both polyline module 1826 and cookie plane modxile 1828. For some 
tj^pes of data, such as inserting a control point, cookie manager 1824 receives data from 
5 polyline module 1826 and passes the data to coolcie plane module 1828. 

Polyline module 1826 manages data related to polyline 1728 and the associated 
control points ui conjunction with polymarker module 1832 and polystate module 1830. 
Polyline 1728 is mainly provided for visual reference. Polyline module 1826 in 
conjunction with polystate module 1830 manages the state of the control points. For 

1 0 instance, in the active state, a control point can be moved or deleted. The active control 

point, such as active control pomt 1720 is preferably highlighted. The control points may 
be enlarged or decreased in size for easier viewing. Polymarker module 1832 provides 
visual context such as highlighting or varied colors, for the control points so that an 
operator Icnows which point is in the active state for moving, deleting, and otherwise 

1 5 editing. Polymarker module 1832 also may provide text such as the location indication 

adjacent active control point 1720. 

Cookie plane module 1828 provides textured geometry, which may for instance be 

lithography-related for geological data, to the surface of the coolcie planes such as cookie 

planes 1712 and 1714. Coolcie state module 1834 monitors die state of the coolcie planes 
20 so that in the active state one or more cookie planes can be moved or otherwise edited 

whereas m the inactive state no changes are made. 

During operation of system 1810, probe module 1822 may notify cookie manager 

1824 that an event has taken place, e.g., a marker deletion, i.e., a control point deletion. 

Cookie manager 1824 then notifies polyline module 1826 which deletes the marker or 
25 control point and joins up the two sxirrounding lines into one line and notifies the 

polymarker module to remove the deleted marker from the list of polymarkers maintained 

by polymarker module 1832. Cookie manager 1824 and cookie plane module then 

convert the two planes into one. 

FIG. 19 illustrates one embodiment of the invention for creating a three 
30 dimensional surface representative of a physical phenomena described by a 3D voliome 

data set such as, for example, a geological fault sxarface. In fact, this method has been 
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found to easily and quickly provide 3D volume visualization for rapid identification and 
interpretation of geological fault surfaces. However, other types of surfaces for other 
types of data could also be quickly described from use of the method of the present 
invention. Although the method is described in more detail below, in general the metliod 
5 is used to digitize control points, such as control points 1902, 1904, 1906, and 1908 along 

a probe face 1910. Visual examination of the textured surface of probe face 1910 permits 
an operator to visually locate control points on a structure of interest such as a suspected 
fault luae. Probe face 1910 is then moved, and a new set of control points are digitized. 
Smface 1912 is then interpolated between the initial points and the new control points. 

10 The control points may be easily edited or moved to more accurately define tlie smface 

may be 1912 at a re£il time frame rate as discussed above. This process can be repeated 
tintil the surface inteipretation is complete, at which time the surface 1912 can be saved. 

Thus, one embodiment of the invention disclosed by FIG. 19 provides a method 
of rapidly constructing a three-dimensional surface or fault found within a 3-D volmne 

15 data set. In a preferred embodiment, 3-D probes, as discussed above, are used. The 

method involves constructing a plurality of spline curves, such as spline curve 1914 on 
the probe face 1910 which may be interactively constructed when visualizing the data 
displayed on probe face 1910. Spline curve 1914 is interpolated using an algorithm 
created in the same way as spline curv^es 1916, 1918, and 1920 are mterpolated. Spline 

20 curves and v-curves 1922, 1924, and 1926, are produced to form a grid. The grid outlines 

the three dimensional surface 1912. 

For construction of the initial spline curve, such as for instance spline curve 1914, 
the user digitizes conti'ol points, such as control points 1902, 1904, 1906 and 1908 on the 
probe face 1910. Markers are produced at these control points and spline curve 1914 is 

25 interpolated between the control points 1902, 1904, 1906 and 1908. Control points 1902, 

1904, 1906 and 1908 may be moved within probe face 1910, thereby interpolating a new 
spline curve 1914. Once the probe face 1910 is selected, the other probe surfaces are 
made transparent for ease of operation. Additionally, the selected probe face 1910 may 
be made opaque in order to view surface 1912 through the probe face 1910. 

30 The user then moves the probe face to 1910 to another position and selects new 

control points. Tlie user may easily move back and forth between previously created 
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5 



spline cxirves by selecting grid intersections such as intersection 1928 or 1930. As 
additional spline curves are created, the v-cur\''es may also be smoothly and quickly 
inteipolated using another algorithm. The user preferably creates a plurality of spline 
curves in the same manner, and interpolation of surface 1912 is immediately displayed 
at real time frame rates as discussed here above. 

The user may stop the probe and move one or more control points, such as control 
points 1902, 1904, 1906 and 1908 to adjust the position of the respective spline curve 
such as splme curve 1914. All other spline curve, such as spline curves 1916, 1918, and 
1920 remain the same while surface 1912 is smoothly interpolated betv^een the current 
10 splinecurve 1914on theprobe face 1910andthepriorsplinecurve 1920. Theremainder 

of the surface 1912 remains the same xmless the user moves the prove face 1910 to 
another splhae curve such as 1920 and proceeds to edit the same thereby reshaping the 
surface 1912 between the current spline curve 1920 on the probe face place 1910 and the 
prior spline cxirve 1918. 

1 5 Additional spline curves may be added between existing spline curves if desired. 

By selecting grid intersections such as grid intersections 1928 or 1930, or by selecting 
control points on the probe face, such as conti'bl points 1902, 1904, 1096 and 1908, the 
user may quickly move the probe face 1910 back and forth as desired. Once a grid 
intersection is selected, the user may move the respective control points within the probe 

20 face 1910 that will be displayed as indicated in FIG. 19, and surface 1912 with its spline 

curves and v-curves will follow interactively. The spluae curves and the v-curves may 
or may not be displayed depending on preference of the user. Preferably, only one probe 
face 1910 is displayed at a time for clarity. 

FIG. 20 discloses modular system 2000 which describes the software for 

25 performing the functions as described in connection with FIG. 1 9. In one embodiment, 

modular system 2000 incorporates many of tlie same modules as used in modular system 
1810 for producing a ribbon section. Thus, the design may permit amore general control 
that may effectively perform both ribbon section functions as well as s\uf ace mapping 
functions. For instance, cookie module 2018, polymarker module 2012, and polyline 

30 module 2014 may be used in the manner discussed above relating to producing ribbon 

sections. Likev^se, probe module 2010 therefore again performs the functions as 
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previously noted such as providing initial sketching splines such as spline 1914, which 
will also be observed to be similar to the polylines discussed hereinbefore. Probe module 
2010 supplies user activity data such as mouse cUcks and keystrokes to tlie various otlier 
modules. Cubic spline module 2016 relates to controlling functions for the editable 
suiface or mapping feature such as surface 1912 shown in FIG. 19. Cubic spline module 
2016, spline manager module 2020, spline curve module 2022, and spline surface module 
2024 are functions that are different from those used only in creating ribbon sections. 
Splme manager 2020 performs a nmnber of different functions such as creating or 
deleting surfaces, changing mode functions from creating a surface to editing a surface, 
as well as read and write functions. For instance^ spline manager 2020 may read and 
write to a surface attribute file regarding attributes such as colors, showing tlie grid, 
marker colors, and the like. Spline curve module 2022 keeps track of the spline curves, 
keeps track of which spline curve may be in the editing or creating state and keeps track 
of changes made to the splme curves. Spline surface module 2024 acts on the changes 
made for interpolating the resulting surface changes. 

By using the system and method of the present invention, geologists and 
geophysicists can more quickly and accurately visualize and interpret 3D seismic data. 
Tliis sharply reduces 3D seismic project cycle time, boosts production from existing 
fields, and finds more reserves. 

Conclusion 

While various embodiments of the present invention have been described above, 
it should be understood tliat they have been presented by way of example only, and not 
limitation. Thus, the breadth and scope of the present invention should not be limited by 
any of the above-described exemplary embodiments, and should be defined in accordance 
Mdfh the spirit of the invention including the following claims and their equivalents. 
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